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ABSTRACT 

An Iterative finite element-integral technique is 
used to predict the sound fiela radiated from the JTlSD 
turbofan inlet. The sound field is divided into two 
regions: the sound field within and near the inlet 
which is c^uteJ using the finite element method and 
the radiation field beyond the inlet which is 
calculated using an integral solution technique. The 
velocity potential formulation of the aco istic wave 
equation was employed in the program. For some single 
mode JTlSD data, the theory and experiment are in good 
agreement for the far field radiation pattern as well 
as suppressor attenuation. Also, the computer program 
is used to simulate flight effects that cannot be 
performed on a ground static test stand. 

NOMENCLATURE 

A area 

^1,2,.. flow coefficients 

C velocity of sound normalized to the stagnation 

speed of sound Cq 

Co stagnation speed of sound 

dg internal duct diameter (normalized by dg to a 

dimensionless value of unity) 

E total number of elements 

e element number 

F function equation, see Eq. (4) 

f initial condition function, see Eq. (14) 

g surface potential in far field, Eq. (16) 

, fr 

k wa.j nuiiiDti' 

N Lota I ' jifiber of unknowns 


Mg Mach number 

m spinning mode number, see Eq. (3) 

N finite element weighting function 

n node number 

p pressure normalized by pgC^ 

R far field position (distance from exit plane 
centerline) normalized by dg 

r radial coordinate normalized by duct diameter dg 
*"o 

t time normalized by dg/Cg 

Z impedance 

exit impedance internal region 
Zg exit impedance external region 

z axial coordinate 

■r specific heat ratio 

e angular coordinate 

p density normalized by pg 

pQ stagnation density 

♦ general flow potential normalized by Cgdg 

a acoustic flow potential 

♦ mean flow potential 

u angular frequency normalized by Cg/dg 

Superscripts 

mean value 


e 


element number 
tfjnjformed value 
SiAscr ipts 

l.j.k.l Index number 
I Iteration Index 

n normal 

s surface 

INTRODUCTION 

For over a decade, the generation, suppression, and 
propagation of noise In turbofan engines has been 
Investigated at the NASA Lewis Research Center. From 
the standpoint of design and capital costs, as well as 
operational manpower, test facilities are expensive to 
operate, lo help reduce future costs, both In-house 
and contractual efforts have been used to develop 
finite difference and finite element programs to 
simulate sound propagation In turbofan engines. 
References (^) and (£) contain a literature review of 
the techniques, advantages, and limitations associated 
with the various msaerlcal solutions of the sound 
propagation equations In ducts. The purpose of the 
present paper is to investigate how accurately a newly 
developed hybrid finite element- Integral program (3) 
simulates the measured noise radiated from the JT15D 
Inlet. 

The velocity potential formulation of the acoustic 
wave equation was employed In the finite 
element-integral program presented herein, since It 
reduced the computer storage and running time by an 
order of magnitude compared to the more general 
linearized gas dynamics formulation. On the other 
hand, when rotational flow exists In the Inlet (wall 
and center body boundary layers), the potential flow 
cannot be used to estimate the effects of shear (4). 
However, since the flow Into an Inlet Is reasonably 
modeled by potential flow and since boundary layers are 
small, the acoustic velocity potential Is Ideally 
suited for acoustic Inlet calculations. 

The turbomachinery noise produced within the engine 
propagates down the inlet duct (also the exit) and 
radiates to the surrounding environment. One major 
problem In numerical simulation of the Inlet Is what 
boundary condition to use at the exit of the duct (^, 
p. 273, or 5). To overcome this problem, the present 
analysis separates the acoustic field Into two regions: 
the "Interior* region within the Inlet and in the vicin- 
ity of the Inlet lip, and the exterior region consist- 
ing of the unbounded region surrounding the Inlet. The 
Interior region contains multidimensional flow while 
the exterior region is assumed to contain only uniform 
flow. In this hybrid method, a finite element program 
Is used to solve for the Interior acoustic field, while 
an Integral technique (b) Is employed In the far 
field. Numerical IteraTlon between the Interior and 
exterior regions Is required to obtain a continuous 
acoustic field across the Interface. 

for comparison between experiment and theory, the 
Inlet data from acoustic treatment experiments with a 
JT150 engine were chosen. In the special JTlbO 
experiment, documented In Ref. (7^), one acoustic mode 
was made to dominate the noise spectrum. Thus, the 
usual problems of noise floors and separation of 
acoustic modes were to a large extent eliminated In the 
data reduction. The results of the present paper will 
Include the effects of the center body. Inlet lip 


shape, and flow gradients on the far field radiation 
pattern and the attenuation of suppressors used In the 
JTISO experiments of Ref. (7). 

Finally, for a static Mgine or fan test rig, one 
major limitation 1$ the Inability to simulate the 
effects of flight on the flow characteristics In the 
Inlet and one the fan source noise and Its radiation. 

A secondary purpose of this paper was to determine the 
changes In che acoustic radiation due to flow entering 
the Inlet as seen by an engine In flight. 

STATIC TEST DATA 

As described in the Introduction, theoretical 
comparisons will be made with the previously published 
Inlet acoustic data for the JT15D engine {Jj. The 
acoustic tests outlined In Ref. (7) were performed at 
the NASA Lewis Vertical Lift FaciTity. A photograph of 
the JTISO engine mounted on the test stand Is shown In 
Fig. 1. The spherical Inflow control screen covering 
the Inlet was used to reduce turbulence and distortion 
In the Inlet flow. The engine exhaust was connected to 
a large muffler to suppress aft fan and Jet noise. The 
engine was mounted 2.9 m or S.4 fan diameters above the 
concrete pad to minimize ground plane effects on the 
Inlet flow. 

The JTISD engine Is a two-spool turbofan engine 
with a nominal 3.3 bypass ratio and a rated thrust of 
9790 N. The fan Is S3.S cm In diameter and has 28 
blades with o6 exit guide vanes. The blade to vane 
number ratio for this fan results In a cutoff 
fundamental rotor-stator Interaction tone. The Inlet 
used with this engine Is shown In Fig. 2. This Inlet 
has a constant 53.3 cm Internal dimneter and Is fitted 
to a bellmouth which In turn Is attached to a large 
constant diameter nacelle. 

A spool piece with 41 equally spaced radial rods 
was attached to the engine front flange, as shown In 
Fig. 2. The wakes of these rods Interact with the 28 
fan blades to produce a BPF tone with 13 
circumferential lobes. At fan speeds between 6400 and 
8400 rpm only the lowest order radial mode can 
propagate, which Is Ideal for the theoretical 
comparisons. 

Provisions were made for an SOOF (single degree of 
freedom) acoustic liner to be placed just upstream of 
the rod section of the Inlet (see Fig. 2). This 
suppressor section was replaceable so that different 
treatment designs could be tested. The acoustic liners 
were designed by theory to attenuate the 13 lobe mode 
produced by the Interaction of the Inlet rods and the 
engine fan. Maximum attenuation efficiency was sought 
for this mode at the blade passage frequency of 3150 Hz 
with an average duct Mach number of -0.147. 

The experimental data will later be presented along 
with the theoretical results. 

MATHEMATICAL FORMULATION 

The formulation of the problem first Involves 
specification of the governing equations and boundary 
conditions. Next, finite element theory will be used 
t, solve for the acoustic parameters within the Inlet, 
while an Integral technique will be employed to 
determine the acoustic parameters outside the engine. 
Two different mathematical techniques are employed here 
to overcome the classical closure problem of specifying 
a nonref lecti ve exit boundary far from the engine. 

Since a direct extension of the finite element 
technique Into the far field Is not possible, the 
integral technique was used In conjunction with the 
finite element algorithm. Iteration will be required 
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to obtain a continuous solution for the acoustic 
quantities at the interface separating the interior and 
exterior domains. Each of these topics will now be 
briefly discussed. 


- -(y - 1) ♦rrV " 2 ^2*2 * ♦r^'" 


-{y - 1) - i [2iii4p] 


( 8 ) 


Acoustic Equations 

Consider the Inlet shown in Fig. 3. The flow of 
air is assumed to be inviscid, nonheat conducting, and 
irrotational. Under these assumptions, the flow in the 
duct can be described in cylindrical coordinates by the 
flowing nondimensional nonlinear partial differential 
equation for the flow potential • (8): 


{Vr * •r''- * * 'll} - •« ■ ^Vu * 

2 2 4 2 

* 2*2*2t * *rVr * * Vzz 

* 

* 

with the Z axis coinciding with the turbofan inlet axis 
and 


C - 1. - (y - 1) 


[*t *7 (*r * ^2) 


where the subscripts indicate partial differentiation 
with respect to the subscripted variables. These 
symbols and those following are listed in the 
nomenclature. 

To obtain the acoustic equations for the inlet, 
first the flow potential a is rewritten a^ the sum of 
a steady axisymmetric mean flow potential ^(r,z) and 
an acoustic potential a 


♦ = a + air.z) 


^-i(iiit - me) 


(3) 


A5 - -{y + 1) - 2 a^^V - (y - 1) a^^a^ 

- (y - 1 ) a^a^/r - i [2^2] ( 9 ) 


^6 “ ~ - w(y - 1)^ Vr * * *zz^ 

where the A's contain all the mean flow quantities. 
The acoustic potential a is complex. 

Since the steady flow is assumed to be 
axisymmetric, the equation expressing conservation of 
acoustic momentum in the e direction can be 
integrated to give the following relationship between 
the pressure perturbation P and the velocity 
potential a: 


p - -p (1*4 + a^a^ a^a) (il) 

Boundary Conditions 

To prescribe the boundary conditions, the engine is 
divided into a finite interior region and an exterior 
region, as shown in Fig. 3 The rear portion of the 
inlet is replaced by an ellipsoidol termination since 
it has been previously established that the geometry of 
this region has very little effect on the sound 
radiated from the inlet (6). 

Along the axis of symmetry, and the acoustically 
hard walls (centerbody, part of the duct wall and the 
nacelle exterior wall) the normal gradient of the 
acoustic velocity potential is 


To account for spinning modes, the acoustic potential 
in £q. (3) has been modified to include e^"*®. 

Also, the solutions have been assumed to be harmonic 
(e~’“t) in time. _ 

The mean flow quantities, ♦, are independent 
of time and thereby satisfy the steady mean flow 
equation. Substituting £q. (3) into Eq. (1), 
eliminating the mean flow quantities, and neglecting 
nonlinear acoustic quantities yields: 


♦n = 0 (12) 

The boundary condition at the surface of a locally 
reacting sound absorbent soft-wall duct can be 
expressed in terms of a specific acoustic impedance. 

Z. Majjigi (9) has shown that a proper formulation of 
the requirement of continuity of particle displacement 
yields 


F(4) = l^i^rr * 


zz 


'^3®rz " 


%®r ^ S®z ^ 




0 (4) 



(13) 


where 





where p is given in Eq. (11). In this analysis, it 
is assumed that the last term in Eq. (13) is negligible 
since ^nn quite small and cannot be a':curaiely 
( 5 \ calculated by the method currently used for computing 
the steady flow. 

The fan plane is the interface between the inlet 
and the turbofan engine noise source (rotor and stator 
blades). Generally, the noise source can be a velocity 
^ ' potential or a pressure source. In the problem 

considered herein a sound pressure source of the form 

( 7 ) P (r) - f (r) (14) 

will be used. As will be shown, a plane radial profile 
(f(r) . constant) and a Bessel function representation 
for f(r) will be used. The pressure P(r) is 
related to *(r) through Eq. (11). 
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Finally, at the Interface between the interior and 
exterior regions an impedance condition 


with a only at the centroid. Thus the total number 
of unknowns M can be expressed as 


^11 “ ^-0 


(15) 


is used to initially close the interior problem. Since 
the value of Zn Is not known a priori, an iteration 
scheme is used to determine the correct value of 
Zii. The subscript I represents an iteration 
counter, while 1 stands for interior region. I 
equals 1 in the first iteration as shown in £q. (15). 

On axis of symnetry in the exterior field, Eq. (12) 
holds again. Along the interface separating the 
interior and exterior regions, the potential in the 
exterior region g(s) must equal the potential in the 
interior region. 

9j(*) “ ♦i(s) (16) 


and the impedances at the interface must be equal 

Zle - Zli (17) 

Again, and gi(s) as well as the impedances 
are not known a priori. The iteration procedure to 
determine these values will be discussed shortly. 

Finally, only outward propagating waves can exist 
in the far field as the distance from the inlet 
approaches infinity. By using the free space Green 
function in the integral solution (6), the absence of 
reflections in the far field is automatically achieved. 

Interior Solution 

A finite element method is used to obtain the 
solutions to Eq. (4) in the interior portion of the 
inlet. To accomplish this task on the JT15D inlet, the 
inlet nas been divided into a number of triangular 
Hermitian elements, as seen in rig. 4. For clarity, 
the number of elements shown have been greatly reduced 
from the actual number used. As shown in rig. 3, the 
indices i,j, and k represent the vertices of the 
triangle and the index i denotes the centroid 
node. The centroid node is not shown in Fig. 4. 

To represent the variation inside the elements of 
the field variable and its derivatives, the 
interpolation function for tne Hermitian elenents is 
written as follows 


♦(r,z) . ^ ^ N^3a,^ - hja^ - N®a,j 

" * '^^zk " "^rk " ^lo*t 


(18) 


where r and z subscripts represent derivatives 
while the i,j,k,i subscripts represent the node 
positions. Values of the weighting functions 

NpN^,... are given in Ref. (10). 

The unknown nodal values of the velocity potential 
a and its derivatives are now constrained to satisfy 
the wave Eq. (4) and the previously discussed boundary 
conditions by using the Galerkin form of tne method of 
weighted residuals. For the simpler Lagrange 
interpolation functions (11) the total number of 
unknowns would be the total number of vertices and 
centroids representing the unknown a values at each 
node. However, using the Hermitian formulation, 
a, ar, and az are unknowns at each vertex 


M - 3x (Number of vertices) E (19) 

where E is the total number of elements each having 
one centroidal unknown. Obviously, the bandwidth of 
the solution matrix increases with the Hermitian 
formulation; however, fewer elements are needed to 
resolve the pressure field. Also, use of the Hermitian 
formulation tends to reduce spurious oscillation 
excited by inaccurate matching at boundaries (^). 

Application of the Galerkin finite element method 
yields the following relationship 

h ff ^ F(4)<1A^*^ - 0 (20) 

e ■ 1 •'•'A 


where F(a) represents the wave Eq. (4) and the area 
of integration is over the element labled e. As given 
in Eq. (18), is the known weighting function 
for each unknown m associated with a node m of 
element e. It should be noted that N^ is zero 
for all elements not having the unknown m (4 or its 
derivative) associated with a particular element. For 
example, vertix node 16 in Fig. 4 will require a 
summation over only six adjacent elements, while the 
vertix node 18 along the boundary will require a 
summation over only three elements. The centroid nodes 
will require only one integration over only a single 
element. 

When Hermitian elements are used, the analytical 
integrations of the element equations required by Eq. 
(20) become formidable. In this case, it is far more 
practical to evaluate the integrals numerically using a 
Gauss quadrature scheme. Since the weighting functions 
are polynomials, the Gauss quadrature is exact (13) and 
the chance of error is substantially reduced by Tetting 
the computer perform the integration. 

Equation (20) provided M equations for the M 
unknown nodal values of 4 and its derivatives. The 
remaining task of assembling and solving the matrix 
equation follows the general approach of finite element 
methods (3 and 10). In this case, however, the 
acoustic velociTTes represented by 4 ^ and 4 j 
are contained in the solution vector and need not be 
calculated by differentiating the potential. 

Exterior Solution 

The integral technique is used to solve the wave 
equation in the far field. In the integral technique, 
the Green's theorem is used to transform the partial 
differential equation within a region to an equivalent 
problem of solving an integral equation over the 
boundary of the region. For acoustic problems, the 
Integrals vanish at infinity and on the axis of 
symmetry leaving only the body surface, or in this case 
the "interface" shown in Fig. 3. The integral 
technique developed by Meyer et al. (6) will be 
employed here. It uses a Green's function that 
satisfies the Somnierf ield condition of no reflection at 
infinity. 

The integral technique can only be applied in cases 
where the properties of the medium are constant or the 
differential equation describing the medium can be 
transformed into a constant property problem. 
Consequently, the interface separating regions I and II 
in Fig. 3 must be extended sufficiently far from the 
exit so that the flow is approximately uniforr, in 
this case, the wave Eq. (4) reouces to 
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2 \ For the calculations reported In this paper, the 

AT \ finite el^nt computer code subdivided the interior 

^ */ region Into 320 Hermitlan elements. These elements use 

’ 1003 nodal points of which 683 are triangle vertices 

(21) with 2049 degrees of freedom (e, v, Ox 

unknowns) and 320 centroid unknowns for a total of 
2369 unknowns. The computer code for the Integral 
technique subdivided the Interface Into 59 segments. 


where 


k - w/c (22) 

and 

Mq - ij/c (23) 

Using the transformations employed by Tester (14) 

r* - r, 2 * - z /yi~ k* - k/^TMT" 

♦ exp ^IkMz*/^ - "o^ ) (2*) 

Eq. (21) becomes 

♦r*r* ^ V/r* ^ * 2 * 2 * " ^ ® 

In Ref. (6), Meyer et al. have transformed t . (2b) 
along with the appropriate boundary conditions Into the 
Integral form. The Meyer et al. approach was used here 
since their computer codes were made available and the 
accuracy of the codes has been previously demonstrated. 

Iteration Scheme 

As mentioned earlier, at the Interface between the 
Interior and exterior regions, the Initial Impedance 
condition, Eq. (15), Is only an approximation. To 
determine the exact value of the exit Impedance as well 
as the acoustic velocities and pressures, a simple 
iterative procedure based on the method of successive 
substitution is formulated as follows: 

1. Employ the exit impedance condition (Z]^^ > 1.) 

and use the finite element method to obtain a solution 
within the interior region. From this solution compute 
the acoustic potential »(s) at the interface. 

2. Let the acoustic potential at the interface as 
computed instep 1 equal g(s). Use the integral 
technique to compute the exterior field. From the 
exterior solution, con^)ute the impedance at the 
interface (i.e.,Zie) where the subscript e denotes 
the external impedance at the Interface. 

3. If Zie is not equal to Zii, set Z^i 
equal to Zi^, and repeat the solution of the 
interior region using the finite element method and 
obtain a new acoustic potential a(s) at the 
interface. 

4. Repeat Steps (2) and (3) until Zji and Zje 
agree within a specified tolerance. 

DISCUSSION OF RESULTS 

As previously outlined, the finite element-integral 
program is now applied to the JTlbD engine 
configuration for the m - 13 mode to determine the 
far field radiation pattern, suppressor effectiveness, 
and effect of flight Mach number. Each of these topics 
will now be considered in order. 


Far Field Radiation Pattern 

The calculated far /(eld radiation pattern for the 
(13,0) mode at 3150 Hz Is shown In Fig. 5 at 24.4 m and 
5.42 m. The similarity between the two far field 
distributions Indlcates^that near field effects are 
absent. The Jump at 30* Is not real and results from 
Insufficient elements. Fewer elements produce a larger 
jump. These results are replotted In Fig. b along with 
the experimental data from Kef. { 7 ), for wnich the 
single m - 13 mode dominates the noise source. The 
data, presented In graphical form In Ref. { 7 )^ are 
tabulated here In Tables I and II for convenience. For 
the noise source at the fan plane, both tne radial 
pressure distribution of the (13,0) mode and also a 
uniform radial pressure distribution In the annulus 
were used. Since only the (13,0) mode propagates, the 
far field pressure distribution was Identical for both 
source Inputs. This, of course, will not be the case 
when more than one spinning mode can propagate. As 
seen In Fig, 6, the analysis and experment agree 
reasonably well. The peaks, however, ere about 10 
degrees apart. 

Also shown in Fig. 6 Is the far field distribution 
for a sharp edged Inlet based on Savkar's (15) 
WIener-Hopf analytical solution. Savkar's THeory was 
modifiea for convection effects using a velocity 
correction factor from Ref. (^) (Eq. (8)). As seen In 
Fig. 6, the sharp lip solution predicts slightly more 
radiation to the higher angles than the thick bellmouth 
finite element solution. 

Soft Kali Ducts 

As mentioned earlier, soft wall suppressors were 
run in the 0T15O inlet to attenuate the sound in the 
far field. The axial position of the suppressor Is 
shown in Fig. 2. In t ‘'2 experiment only half of the 
suppressor was active. As seen by the slanted lines 
originating at nodes 18 and 44 in Fig. 4, tne finite 
element nodes were positioned exactly at the entrance 
and exit of the active liner. The second half of the 
suppressor was taped to give a hard wall. With the 
full suppressor open, the tone level fell below the 
broadband noise level which prevented a quantitative 
evaluation of the suppressor's performance. Three 
suppressor designs were built having estimated nominal 
values of resistance of 2.272, 1.136, and 0.638 (Taole 
1, (7^)). The suppressors were all specially built to 
have identical reactances of 0.5 which is stiffness 
controlled at 3150 Hz. 

, Figure 7 displays the calculated far field 
radiation pattern for the same source as in Fig. 5 but 
with a soft wall suppressor having a resistance of 
0.638 and reactance of 0.5. A comparison of Fig. 5 
with Fig. 7 shows the sound has the same far field 
shape but is decreased in magnitude by about 22 88. 

A comparison with the experimental data will now be 
made with data between 50 and 80 degrees. These angles 
were chosen since in the experiment the sound pressure 
level for these angles stood well above the broadband 
noise floor. At 60 degrees, subtracting the SPL of 
69.02 in Fig. 7 from the haro wall value of 91.23 In 
Fig. 5 yields a calculated attenuation of 22.21 d8. 

The other angles also have the same theoretical 
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attenuation. This attenuation along with other 
calculated attenuations for the other resistance values 
are plotted in Fig. 8. As seen In Fig. 8, the 
xperiment and theory are again in reasonable agreement. 

Also bhown in Fig. 8 is the attenuation calculated 
using a semi- infinite soft wall duct theory (solid 
circular point) which considers only the lowest radial 
mode for m > 13. The maximum attenuation for the semi- 
infinite duct analysis is 30 dB ((2), pg. 8 or Fig. 9) 
which is approximately 6 dB higher” than the peak 
numerical prediction and about 13 dB higher than the 
numerical prediction at a - 1.14. 

Note also that the numerically determined maximum 
attenuation (at x - 0.5) occurs at a lower resistance 
(^0.5) than that determined using a single radial mode 
(e~1.14). This difference can be explained using 
the argument presented in Ref. (7) tnat modal 
scattering is important for this”extreme case of a 
single source mode entering a very short acoustic liner 
(L/D - 0.15). The impedance discontinuity causes the 
single hardwall mode to scatter into several radial 
modes in the soft wall section. Since the liner is so 
short the higher radial modes contribute significantly 
to the overall suppression. It was shown in Ref. ( 17 ) 
that higher radial modes have progressively lower 
optimum resistance. It would thus be expected that a 
multimodal liner solution (which the numerical solution 
is) would produce a lower optimum resistance than the 
lowest radial mode solution used 1. design these 
liners. This explanation of modal scattering being 
significant is reinforced by the calculations shown in 
Fig. 9. As the acoustic liner is lengthened (smaller 
0/L), the numerically calculated attenuation approaches 
the single mode attenuation shown at 0/L > 0. For an 
extremely long duct this lowest mode would uominate the 
sound attenuation. 

Flight Mach Number 

As mentioned In the introduction, static 
experiments on test stands as shown in Fig. 10(a) 
cannot simulate the true flight situation as shown in 
Fig. 10(b). Although a bellniouth inlet is not a 
realistic inlet to fly, a calculation will now be 
performed to demonstrate the capability of the 
program. Figure 11 shows the far field radiation 
pattern for the case of flight at 44.7 m/sec (100 
miles/hr). A comparison to the static case of Fig. 5 
is shown in Fig. 12. As seen in Fig. 12, the noise is 
shifted by flight to lower angles although the peak did 
not seem to move. These results cannot be generalized 
to flight, because a bellmouth inlet would not be used 
in flight. 

LIMITATIONS OF THEORY 

As seen in the previous section, the theory is 
versatile and gives reasonable predictions of the far 
fielo radiation pattern and suppressor attenuation. 

This data was for a nearly cutoff mode with a 
relatively long wave length. For well propagating 
mooes with relatively short axial wavelength, however, 
the present computer does not have the storage capacity 
to resolve the many axial wavelengths. This is 
particularly true when the Mach approaches unity and 
the convective effects further shorten the sound 
wavelength. Therefore, research is continuing (2) 
using special numerical techniques such as the wave 
envelope concept (jS and 19) to reduce the required 
computer storage and the Transient technique (^) to 
completely eliminate matrix storage requirements. 

Finally, for well propagating modes the boundary 
layer can have significant effects The present 


potential flow formulation does not Include the effects 
of a boundary layer. In this situation, the complete 
solution to the separate linearized continuity, 
momentum, and energy equations will be required. With 
this increased number of unknowns, the wave envelope 
technique or a transient approach (21 and 22) looks 
very attractive. 

CONCLUSIONS 

Finite element theory coupled with the integral 
technique is a versatile tool for handling noise 
propagation in engine inlets and radiation to the far 
field. The theory can also be applied to any situation 
where the axial wavelength of sound is relatively 
large, such as an inlet with low frequency sound 
excitation and low Mach number and an exhause duct 
where convective effects tend to increase the sound 
wavelength. 

For the single mode JT150 data, the theory and 
experiment are in reasonably good agreement for the far 
field pattern as well as the suppressor attenuation. 
Flight effects are also easily handled by the program. 
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TABLE I 

HARD WALL JT15D TEST INLET - Fig. 6 
(Far field 24.38 m, 3150 Hz, m • 13, Ref. (2)) 

SPL 

Angle from inlet, degrees 



40 

50 

60 

70 

80 

90 

Run 1 

88.5 

94.0 

95.6 

95.4 

90.0 

85.5 

Run 2 

87.2 

94.4 

96.4 

96.0 

90.5 

87.7 

Run 3 

88.5 

96.3 

96.0 

95.7 

89.5 

84.6 

Run 4 

88.1 

96.2 

97.7 

96.3 

90.2 

85.0 

Run 5 

88.8 

96.7 

97.8 

96.1 

90.0 

86.4 

Run 6 

88.5 

96.3 

96.7 

95.2 

90.3 

85.7 

Average 

88.3 

95.7 

96.7 

95.8 

90.1 

85.8 

Normalized 

91.6 

99.0 

100.0 

99.1 

93.4 

89.1 


TABLE II 

FAR FIELD SOFT WALL JT150 SPL ATTENUATION - Fig. 8 



(x 

« 0.5 1 R6f • 

0) 

9 ’ 

SPL 


Wall 

Runs 

Angle from 

inlet, degrees 

Resistance 


50 

60 

70 

80 

2.272 

Run 1 

8.5 

9.6 

9.6 

9.0 


Run 2 

7.4 

8.8 

10.7 

10.0 


Average 

8.0 

9.2 

10.2 

9.5 

1.136 

Run 3 

19.2 

17.6 

20.7 

18.3 


Run 4 

16.2 

17.9 

20.3 

20.8 


Average 

17.7 

17.8 

20.5 

19.6 

0.638 

Run 5 

23.7 

24.0 

21.6 

19.0 


Run 6 

22.2 

20.6 

20.0 

20.5 


Average 

23.0 

22.3 

20.8 

19.8 
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Figure 1. - JT15D engine with inflow control device on test stand with exhaust muffler. 




'“INFLOW CONTROL DEVICE 
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Figure Z - Test configuration of JT 15 D engine inlet used for single mode SDOF suppressor experiments, dimensions in inches (cm). 
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Figure 3. - Geometry of an idealized inlet 
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Figure 5. - Theoretical far field SPL directionality pattern 
for 3150 Hz and (13, « mode at source for hard wall con- 
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(D) Flight test 

Figure 10. - Effect of forward velocity on iniet stream lines. 
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Figure IL - Theoretical far field SPL directionality pat- 
tern for flight at 44, 7 m/sec OOO mlles/hrt for 3150 Hz 
and (13, 01 plane mode at source for hard wall config- 
uration. 



